


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1980-09 


Heat exchanger optimization. 


Hedderich, Conrad P. 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/17604 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
i (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


HEAT EXCHANGER OPTIMIZATION 


Conrad P, Hedderich 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





PeaeslS 


HEAT EXCHANGER OPTIMIZATION 


by 
Conrad P. Hedderich 


September 1980 


iThesis Co-Advisors: M. Kelleher 
G. Vanderolaats 





Approved for public release;distribution unlimited. 


7197430 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


. REPORT NUMBER 


4. TITLE (and Subtitie) 
Heat Exchanger Optimization 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


2. GOVT ACCESSION NO, 3. RECIPIENT'S CATALOG NUMBER 












§. TYPE OF REPORT & PERIOD COVEREO 
Master's Thesis; 
September 1980 


6. PERFORMING ORG. REPORT NUMBER 










7. AU THOR(a) . CONTRACT OR GRANT NUMBER(S) 






Conrad P. Hedderich 










- PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 







9. PERFORMING ORGANIZATION NAME ANO AOORESS 


Naval Postgraduate School 
Monterey, California 93940 



















REPORT OATE 


11. CONTROLLING OFFICE NAME ANO ADORESS 


Naval Postgraduate School 
Monterey, California 93940 





188 


. MONITORING AGENCY NAME 4& ADORESS(/f different from Controlling Office) | 15. SECURITY CLASS. (of thie report) — 


Unclassified 


Se. OECLASSIFICATION/ COWNGRADING 
SCHEOULE 


» OISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited. 


. DISTRIBUTION STATEMENT (of the abetraci entered in Block 20, tf different from Report) 


- SUPPLEMENTARY NOTES 


. KEY WORDS (Continue on reverse eide if neceesary and identity dy diock number) 


Heat Exchanger Design Using Numerical Optimization 
HEDSUP Design 


Heat Exchanger Spit rmescieLon 
». ABSTRACT (Continue on reverse side if neceesary and identify by block number) 


A computer code was developed for the analysis of air-cooled 
heat exchangers and was coupled with a constrained function mini- 
mization code to produce an automated air-colled heat exchanger 
design and optimization program with many new capabilities. 








A general iteration free approximation method was used for 
the analysis which calculates the mean overall heat transfer 
coefficient and the overall pressure drop for many flow 






DD , Seuies 1473. —s- EDITION OF | NOV 6318 OBSOLETE UNCLASSIFIED 


S/N 0102-014* 6601 | =e 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 





FOCUINTY CLASEIFIC ATION OF TuIS PAGE/When More Batered. 





#20 - ABSTRACT - (CONTINUED) 


arrangements, taking into account the variation of the heat 
transfer coefficients and the pressure drop with temperature 
and/or length of flow path. 


The code is not limited to surfaces found in the litera- 
ture, but will accommodate any triangular pitch bank of finned 
tubes in multiple-pass configurations. 


a pon, 1463 UNCLASSIFIED 
S/ 0102-014-6601 2 SECURITY CLABBIFIC ATION QF THIS PAGESWrhenm Seaia Eniored) 





Approved for public release; distribution unlimited. 


HEAT EXCHANGER OPTIMIZATION 


by 


Conrader. sledge mich 
Lieutenant, United States Navy 
B.S.M.E., United States Naval Academy, 1973 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 
NAVAL POSTGRADUATE SCHOOL 


September 1980 





ABSTRACT 


A computer code was developed for the analysis of air- 
cooled heat exchangers and was coupled with a constrained 
function minimization code to produce an automated air-cooled 
heat exchanger design and optimization program with many 
new capabilities. 

A general iteration free approximation method was used 
for the analysis which calculates the mean overall heat 
transfer coefficient and the overall pressure drop for many 
flow arrangements, taking into account the variation of the 
heat transfer coefficients and the pressure drop with tempera- 
ture and/or length of flow path. 

The code is not limited to surfaces found in the litera- 
ture, but will accommodate any triangular pitch bank of finned 


tubes in multiple-pass configurations. 
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NOMENCLATURE 


English Letter Symbols 


A - total heat transfer area, i 


e =- specific heat, BTU/lbm-°F 

e - heat capacity rate, BTU/hr-°F = m Cy 

D - diameter, in. 

E euetwaetion factor 

i =e) correction factor 

a peceleration of gravity, 32.2 ft/sec 

H - corrected heat transfer coefficient, BTU/hr-f£t*-°F 
h - bank height, in. 

J mecolourn factor 

k Seecnermal conductivity, BTU/hr-f£t-°F 

L - length, in. 

Q aera height, in. 

m - mass flow rate, lbm/hr 

m - 2H 7k,t , ft 7 

n Seeoumber Of $ (used with appropriate subscript) 
N Seeilmber Of tubes 


p - pressure, psi 
P mel cch, in. 


Pr - Prandtl number 


See heat transfer rate, Btu/hr 
a - radius, in. 
R - heat transfer resistance, hr-£t¢-°F/BTU 





Re 


Reynolds number 


distance between adjacent fins, in. 


fin spacing center-to-center, 


fin thickness, in. 


true temperature, deg. 


overall heat transfer coefficient, 


bank width, in. 


Greek Letter Symbols 


AT - temperature difference, deg. 
n - surface efficiency 

$ - fin efficiency 

u seeeviscosity, lbm/ft—hr 
0 - density, lbm/£t> 

W SecemDeracure correc 
Subscripts 

a ae air 

c Saeco ld 

Hf eee 12 

moe cree face 

h Se @ to 

al = inside 

5 - reference number, 

Q = limiting 

L - tliongitudinal 

m - mean 

O - outside 


E10n 


iE ang Stik 


ttle 


BTU/hr-f£t*- °F 





p - pass 


a - rows 
fe Secrahioverse 

a = total 

W Se wWacer 

X ~ cross~sectional (flow) 
a - entering 

2 - leaving 

co - ambient 


Eueerscripts 


o - uncorrected 
* - estimated 
alt See initial 
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iN Roo UeT LON 


A. BACKGROUND 
The cooling of fluids by passing ambient air over extended 
tube surfaces is a relatively recent development in heat 
exchangers. Its application has come about cautiously, due 
to the usual reluctance to change from well established and 
well documented methods, i.e. the shell-and-tube heat exchanger. 
However, concern for the environment and economic pres- 
Sures have necesSitated the use of air aS a coolant. Smith 
[1] has listed some typical advantages of direct cooling with 
air as compared to cooling with water in a shell-and-tube 
exchanger: 


a. Eliminates the problem of temperature rise in, and 
pollution of, water resources. 


b. Enables plant location to be independent of a water 
supply. : 


c. Eliminates the necessity of much coolant piping. 


d. Reduces heat exchanger maintenance costs by eliminating 


the need of descaling water-side surfaces. The 
mechanical drives will operate in a noncorrosive 
atmosphere. 


e. Eliminates water treatment. 

f. Limiting coolant temperatures is unnecesSary. 

g. Enables installation of exchangers at elevations above 
other operating equipment at no penalty, thus reducing 
ground area requirements. 

With air colling becoming more and more competitive with 


water; even when water supplies are plentiful, an automated 


air-cooled heat exchanger design package (which could be used 
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for trade-off studies, first cut analysis, and conceptual 
design), would be of great use. 

The design of an air-cooled cross flow heat exchanger is 
a complex task requiring the examination and optimization of 
a wide variety of heat transfer surfaces. Studies have shown 
that a poor choice of either the heat transfer surfaces or 
design parameters can more than double the costs charageable 
to a heat exchanger [2]. 

For the optimized design of heat exchangers with the 
computer; reliable, but fast, calculation methods for the 
mean overall heat transfer coefficient and the overall 
pressure drops are needed for the following reasons: 

a. Conventional simple methods using mean values of 
temperatures as reference temperatures can lead to 
undesirable errors [3]. 

I. Numerical stepwise integrations are prohibitively 
time consuming. 

Boe ROV LEW 

A numberof heat exchanger design methods have been proposed 
to determine the optimum heat exchanger design. Bergles, et 
al. [4], performed an evaluation of different objective func- 
tions for compact heat exchangers with different heat trans- 
Fer surfaces, but the same specifications. The method did not 
include any actual optimization techniques, but results did 
Show that a great improvement in heat exchanger performance 


can be made by proper selection of design parameters. 
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The method of Fax and Mills [5], used Lagrange multi- 
pliers to optimize a heat exchanger design under specified 
constraints. This technique required that the objective 
function and constraints be expressed explicitly and be differ- 
entiable throughout the range of interest. The total num- 
ber of constraints had to be less than the total number of 
variables, and all constraints had to be equality constraints. 
The method was obviously restricted to a very limited number 
of problems. 

For a unit section of an air-cooled exchanger of standard 
length and width equipped with specified fans, both Schoonman 
[6], and Joyce [7], used factorial searches of fin spacing, 
number of rows and air rate to maximize the ratio of heat 
transfer to cost. Because the number of sections must be 
rounded to an integer, this method only appears suitable for 
large exchangers. 

Nakayama [3], used a similar approach requiring the plotting 
of heat transfer coefficients and pressure drops which pre- 
sumably could be programmed into a factorial search. 

Kern [8] derived an analytic expression for the annual 
cost as a function of air rate and number of rows. (This was 
done by uSing a constant tubeside heat transfer coefficient 
and an arithmetic mean temperature difference.) The optimum 
was found by setting partial derivatives equal to zero. 
Alternatively, if the airside is assumed to control, the 


optimum allocation of the components of total cost might be 
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found through geometric programming, as illustrated by Auriel 
and Wilde [9]. Oshwald and Kochenberger [10], also presented 
a geometric programming method for heat exchanger optimization 
and used it to select heat exchanger fluids considering power 
requirements, cost, tube diameters, velocities, temperature, 
and other physical properties. 

After discussing factorial, univariate, and random search 
meaenods Cor an optimum design of a shell-and-tube exchanger, 
Briggs and Evans [{11, 12], discuss a "logical search method"; 
or what Peters and Nicole [13], call "heuristics". With this 
method, an engineer makes use of selecting design variables 
close to the optimum to obtain an optimum design. The heuris- 
tic method is less scientific, and is useful only when computer 
time and storage are at a premium. Due to the large number 
of discrete variables encountered in air-cooled heat exchanger 
design, Peters and Nicole [13], chose to base their cost- 
optimizing design programs on heuristic algorithms (starting 


close to optimum), specific to the equipment under considera- 


e1on . 
Mott, et al. [2], discuss a computerized procedure for 
designing a minimum cost heat exchanger. The method minimizes 


a cost index expressed as a function of fluid pumping power. 
The algorithm imposed no constraints. 

To this point, no mention has been made of applying the 
concepts and techniques of nonlinear programming to optimizing 
the design of heat exchangers. However, Palen, et al. [14], 


in 1974 proposed using the Complex Method [15], for the heat 
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exchanger optimization problem. They found a minimum cost 
shell-and-tube exchanger by varying six geometrical parameters. 
The Complex Method requires several feasible starting designs 
before optimization can be performed. 

Johnson, et al. [16], coupled an existing shell-and-tube 
condenser design code with a constrained function minimization 
code to produce an automated marine condenser design program 
of vastly different complexity. 

The most complete work to date has been accomplished by 
Afimiwala [17]. He has applied various nonlinear programming 
methods of optimization to the heat exchanger design problem; 
including an experimental interactive graphical approach and 
exterior penalty function techniques. The gradient based 
search methods of Davidson - Fletcher - Powell and conjugate 
gradient were used for the resulting unconstrained minimiza- 
tions. The exterior penalty method is extremely useful, since 
an initial solution satisfying the constraints is not required. 
The gradient based search methods are efficient when considering 
computer time. 

Finally, Fontein and Wassink [18] utilized the complex 
method of Nelder and Mead [19], and a steepest descent method 
[20], for optimizing a shell-and-tube exchanger. 

It can be seen that although there are many methods that 
have been presented for heat exchanger optimization, each of 
the methods has its own limitations; none is completely general. 
Of all the design procedures cited above (those of which are 


applicable to cross flow air-cooled heat exchangers), all are 
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limited to the 120 individual surfaces found in the open 
literature [21] for the calculation of the air-side heat 
Peanster coefficient and friction factor. Therefore, the 
designer is faced with choosing an optimum surface from a 
number of individual optimal designs calculated from one of 

the above methods. In addition, the above methods treat the 
overall heat transfer coefficient as a constant, or they become 
involved with time-consuming numerical stepwise integrations 

in an attempt to account for the varying heat transfer coeffi- 
cients. 

This paper tries to bridge this gap by presenting an 
optimization routine that: selects an optimal surface, takes 
into account the varying heat transfer coefficients and fric- 
tion factors across the exchanger, performs each analysis in 
an iterative-free manner, and may start with an infeasible 


design. 


C. METHODOLOGY 

faethe Control Program for Engineering Synthesis and 
Boenstrained Function Minimization (COPES/CONMIN) optimizing 
scheme, a nonlinear optimization program is available that is 
capable of optimizing a wide class of engineering problems 
[22,23]. Therefore, for the heat exchanger design problem, 
it was necessary to develop a subroutine, which given a start- 
ing design, would analyze an air-cooled heat exchanger, and 
which would be suitable for coupling with the optimizer, 


COPES/CONMIN. 
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Figure 1 illustrates the procedure by which the heat 
exchanger was analyzed. Initial input consisted of a com- 
plete listing of all design parameters, whether known or 
estimated. Those that were estimates, 1.e. unknowns, were 
later passed to the optimizer as design variables. 

The analysis scheme and the optimizer will be discussed 


in much greater detail in the following chapters. 


Pee OBJECTIVES 

The objectives of this thesis are two-fold. 

The first objective was to develop a computer code, 
hereafter referred to as ANALIZ, which would analyze an 
air-cooled heat exchanger given any initial design. The 
analysis scheme was to: be iterative-free, take into account 
the variation of the heat transfer coefficients and the pres- 
sure drop with temperature and/or length of flow path, and 
finally, be written in such a manner that it could be coupled 
with an existing optimizer COPES/CONMIN. 

The second objective was to actually couple ANALIZ with 
a numerical optimization program. This would produce a 
detailed design program which would have the capability 
to determine an optimum surface, while optimizing the objec- 


tive function such as size, weight, cost, etc. 
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IIT. NUMERICAL OPTIMIZATION 


A. BACKGROUND 

Almost all design problems require either the maximiza- 
tion or the minimization of some parameter or function. 

This parameter shall be called the design's objective function 
[24]. For example, the problem may call for a heat exchanger 
with a minimum volume. The expression for volume would be 

the design objective function. For the design to be accepta- 
ble, it must satisfy certain design constraints. For exam- 
ple, an air heater must be designed so that it will fit into 

a given space. Therefore, the engineer must set design con- 
straints on the maximum size of the exchanger. 

If the objective function could be easily formulated 
analytically, the maxima or minima could be found by using 
the methods of differential calculus. However, the limita- 
tions of this method are obvious. 

Another numerical method that would be satisfactory for 
small scale problems would be an iterative solution technique. 
A computer program could be written containing a series of 
nested iteration loops that would vary the design parameters 
and solve the problem for a variety of values for each of the 
parameters. For other than small, easily formulated problems, 
the cost in central processor (CPU), time would be prohibitive. 

Over the last twenty years, many numerical optimization 


techniques have been developed specifically for computer 
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utilization. These techniques usually do not require a 
specific algebraic equation, but rather any computer algorithm 
to which design variables can be input and from which the 
objective function and design constraint values can be 
determined is acceptable. For this reason, nonlinear pro- 
gramming methods were chosen for the air-cooled heat exchanger 
design. Some of these techniques were summarized by Shah et 
miei. 25) ian figure 2. 


1. One-Dimensional Search Methods. Two of the most common 
of these search methods are the golden section [26] 
and quadratic interpolation [27]. The former iso- 
lates the minimum in regions of successively decreasing 
Size, the latter performs a series of iterations 
approximating the objective function as a quadratic. 


2. Multidimensional Unconstrained Search Methods. These 
unconstrained searches can be performed by a sequence 
of one-dimensional minimizations in the proper direc- 
tions. 


3. Multidimensional Constrained Search Methods. A common 
method here for enforcing the constraints in an optimi- 
zation scheme are based on the sequential penalty func- 
tion method. These techniques convert the constrained 
optimization problem into a sequence of unconstrained 
problems. This is accomplished by applying either an 
exterior or interior penalty to the objective function. 
The Complex Method [15], locates the optimum based on 
an intuitive approach in a n-dimensional space defined 
by the independent design variables. The method of 
feasible directions is used primarily for inequality 
constraints. 


An optimization program based upon the Augmented Lagrangian 
Multiplier Method and the method of feasible directions was 


chosen for this research project. 


B. CONSTRAINED FUNCTION MINIMIZATION (CONMIN) 
Vanderplaats [22], developed an optimization program, 


CONMIN, based on the method of feasible directions which is 
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capable of optimizing a wide variety of engineering problems. 
CONMIN is a FORTRAN program, in subprogram form, that opti- 
mizes a function subject to a set of inequality constraints. 

The following definitions will be useful in the following 
discussion: 


ime besign Variables — those parameters which the optimiza- 
tion program can change in order to improve the design. 


mee Design Constraints - those parameters which must not 
exceed given bounds for the design to be acceptable. 


Pee jceccive Function — the parameter which is going to 
be minimized or maximized. 


The general nonlinear inequality constrained optimization 


problem can be written mathematically as follows [28]: 


Minimize F(X) (1) 


Subject to: 


oe (X) < 0 j = 1,NCON (2) 
a nx k= 1 NDV ay 
ile = a — Ac 
where 
a 
- . 
x xX, (4) 
Aupv 


The vector X is the vector of design variables, with NDV 


equal to the number of design variables. The objective function, 
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F(X), given by eq. (1), as well as the constraint functions 
@ieven by eq. (2), may be linear or nonlinear functions of 
the design variables. They shall also be explicit or implicit 
functions of X, but must have continuous first derivatives. 
NCON is the number of constraints. x; and x are the lower 
and upper bounds or side constraints placed on the design 
Variables. Side constraints could be included in eq. (2), 
but are treated separately for efficiency. Equality con- 
straints are not dealt with by CONMIN, but will be treated 
separately by a multiplier method. 

CONMIN requires that an initial vector of design variables, 


meewhich may or may not yield a feasible design, be specified. 


The design process continues iteratively as: 


xatl 2 XS gy ge SF (5) 


where S72 is a vector search direction, a* 1S a scalar quantity 
which determines the amount of change in X and gq is the itera- 
tion number. At iteration g a direction S7 must be found 

which will reduce the objective (usable direction for mini- 
mization), without violating any constraints (feasible direction), 
see Figure 3 [30]. Once 8S? is determined, eq. (5) becomes a 
one-dimensional search problem in whcih o* must be found 

such that; F(X) is at a minimum (see Figure 4 [29]), a new 
constraint 1s encountered, or a currently active constraint 

(9, (X) = 0) is encountered again. 


The design problem at iteration gq+l becomes one of 


finding a usable-feasible direction, ae and a move parameter 
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a*., This process is illustrated geometrically by Johnson 
[30]. Consider a condenser problem with just two design 


variables, xX, and Xo, where 


Xy = condenser tube outside diameter, 


X tube pitch to diameter ratio. 


2 


Let the objective function be condenser volume, VOL (xX). 
Assume that the tube bundle diameter must be greater than a 


given value, BD win’ and that the cooling water pumping power 


= 


must be less than a given value HP ax’ PiIguse. oS lwiwstrates 
the design problem geometrically. 

It should be reiterated here, that while Johnson's exam- 
ple starts with a feasible initial design, (a) Pe elds mes eine 
a requirement and CONMIN is capable of optimizing given an 
infeasible initial design. This is obviously a great benefit. 

The optimization process begins by calculating the gradient 
Baetene Ob Jective function by using finite difference. Each 
design variable is perturbed by .0l in a single forward step. 

The gradient of the objective function, VF, shown in 
Figure 6 is simply the vector of the first partial deriva- 


tives with respect to the design variables; that is: 





oF AVOL 

a - IX) . AX, 

VF(X) = VVOL = es (6) 
oF AVOL 





om tl pe | 
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Therefore, because no constraints are active or violated at 
(A) . VF defines the direction of steepest ascent. Because 
it is desired to minimize F, the greatest improvement can be 


Made by moving in the negative gradient direction so that 
= VOL. ce 


With the value of S now determined, a search is per- 
formed from (a) until the minimum F is found at on Figure 
6. This is accomplished by taking several values of cl in 
eee) and interpolating for the a* which will give the mini- 
mum value of F. 

The second design iteration is begun at by again 
perturbing X to find VF. Instead of moving in the direction 
of steepest descent, a conjugate direction, developed by 
Fletcher and Reeves [31], is chosen with this method, S is 
Calculated as follows: 


gO 
a | VE (XK) 2 | q-1 
Sa = _- VF (xX) 7 + yx) te = ‘ 


(3) 
see Figure 7. 

The Fletcher-Reeves method is used in order to speed 
convergence. With the new conjugate direction, a search is 
performed in this direction until a constraint is encountered. 
This occurs at (C) of Figure 8 on the pumping power constraint. 

At (Cc) Pov oeheie nC HP constraint active, not only is VF 
found, but the gradient of the active constraint is also 


computed, again using finite difference. The requirements on 
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Bae mew Search direction are now twofold; it must reduce the 
objective function and, at the same time, not violate the 
active constraint. This is solved by using the method of 
feasible directions developed by Zoutendijk [32], and imple- 
mented by Vanderplaats and Moses [33]. 


The problem of finding the new S can be stated as [34]: 


Maximize g subject to the constraints 


Pee Set 6 = 0 (2) 

Vg, (X)°S + Osea < ae j = 1,NAC (0) 

eo cmc (11) 
where Vg, (X) = <- VHP and NAC is the number of active constraints 


(in this case NAC = l). 

If equation (9) is satisfied and g iS positive, the search 
direction will reduce the objective function and is defined 
as a usable. If equation (10) is satisfied and g 1s positive, 
Sis a feasible direction, because in this direction, no con- 
Straints will be violated if only a small move is taken. See 
Figure 3. O. is defined as the push-off factor for the active 
constraint and causes the design to move away from the con- 
Seraint. o. must be greater than or equal to zero in order 
to maintain a feasible design. If the maximum value of 8 from 
eemaetons (9) through (11) is zero, then there is no direction 
that will both reduce the objective function and also be feasi- 


ble. Therefore, the current design is, at least, a local 
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minimum. In Johnson's example, a usable-feasible direction 
exists and a one-dimensional search leads to © in Figure 8 
where the minimum bundle diameter, BD in’ constraint is met. 

From (0) it should be noted here that CONMIN had informa- 
tion regarding the linearity of the BD in Constraine- and, 
therefore, in (10), has set oF = 0 to allow S to follow the 
constraint as shown in Figure 9. The one-dimensional search 
along this constraint is carried out until no further design 
improvement is realized. This occurs at () : 

This discussion of CONMIN would not be complete without 
Citing the program's limitations. NDV directly affects the 
computational time required to reach the optimum. Since the 
calculation of gradients required for each design variable at 
the beginning of each design iteration is found by a finite 
difference step, which requires a complete pass through the 
analysis portion of the program, there is a subsequent in- 
crease in CPU time as NDV increases. Also, due to the inter- 
action between design variables, as NDV gets larger, convergence 
slows during the optimization process. Vanderplaats [24], 
recommends that for most problems of general interest a prac- 
tical limit of NDV = 20 be imposed. NCON does not present 
the same problem because gradient information is calculated 
Simultaneously with VF and then only if the constraint is 
active or violated. 

CONMIN offers no guarantees that a global minimum has been 


reached. Therefore, to lend some assurance, the design is 
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started with several different initial vectors until the 
same optimal design is reached. 
Although CONMIN performs very well with inequality con- 


straints, equality constraints such as: 


cannot be dealt with directly, but must be treated separately, 
using a different method, which will be discussed in the 


following sections. 


C. CONTROL PROGRAM FOR ENGINEERING SYNTHESIS (COPES) 

Recall that CONMIN was written in subroutine form, 
Vanderplaats [23], has developed a main program which greatly 
enhances the use of CONMIN. 

For this main program, COPES, the user must supply an 
analysis subroutine titled ANALIZ. The subroutines CONMIN 
and ANALIZ are then used by COPES to optimize the objective 
mimetion subject only to the inequality constraints. 

ANALIZ must be organized into three segments: input, 
analysis and output. Based on the value of a counter, ICALC, 
ANALIZ performs the proper function in sequence. 

The COPES program currently provides four specific capa- 
bilities; 


1. Single Analysis - one cycle through the program, as 
if ANALIZ was executing alone. 


2. Optimization - minimization or maximization of the 


objective function with constraints and side- 
constraints imposed. 
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3. Sensitivity Analysis - used to explore the effect of 
changing one or more design variables on one or more 
FUNG Eons. 


4. Two Variable Function Space - provides tables of data 
of all specified combination of two design variables. 


However, a recent addition to COPES (still in the develop- 
femeal Stages), has put the Augmented Lagrangian Multiplier 
Method (ALMM) at the disposal of the programmer. Because of 
its good rate of convergence and its theoretical properties, 
the ALMM is preferred for equality constrained problems [35]. 

Therefore, COPES will call the various subroutines in order 
to optimize the objective function; subroutine ANALIZ for 
necessary analysis information, CONMIN for an optimum based 
On inequality constraints, and a subroutine (yet to be named), 
utilizing ALMM for an optimum satisfying the equality con- 
straints. 

For detailed explanation of the ALMM, its background and 
mathematical derivation, consult reference [35]; but for now 


consider the equality constrained problem: 
Min £(X) (10a) 


Spubgectetoe 94 (X) =| © k = 1,NECON (10b) 


| 
where NECON is the number of equality constraints. 
Define the "modified" Lagrangian function as: 
— = — m — 
L(,%) =  £(X) + a A. hy (X) 


where hs is the Lagrangian multiplier. The problem can now 


be stated as: 
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Min Tex) (10c) 


Subject to hy (X) = 0 (10d) 


Then, according to Lagrange, if a \ can be found for which 
X solves the problem stated above, then X is also the solution 
to the original problem, eq. (10a) and (10b). 

The new problem is solved by the conventional exterior 
penalty function method, because this is believed to be one 
Snetne most efficient algorithms for the solution of such 


equality constrained problems [35]. 
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II. HEAT EXCHANGER ANALYSIS 


A. INTRODUCTION 
To meet the objectives of this thesis, an analysis program 
for an air-cooled heat exchanger must be coupled with a 
numerical optimization scheme to produce a complete, detailed 
design package. COPES/CONMIN has greatly simplified this task. 
This analysis program must be written in subroutine form, 
titled ANALIZ, and organized into three segments: input, 
execution, and output. The analysis subroutine must also: 
Mmeltake into account the variation of the heat-transfer 
coefficients and differential pressure drop with 
temperature and/or length of flow path. 


2. Be iterative free, if possible. 


3. Be written in such a manner that the optimizer will 
play a role in surface selection. 


With the number of design variables approaching the prac- 
tical limit, the importance of an iterative free analysis 
subroutine cannot be over-emphasized. The reason being, that 
at the beginning of each design iteration in CONMIN, the calcu- 
lation of all gradients (each design variable and active con- 
straint) requires a complete pass through ANALIZ. Therefore, 
the computational time required by ANALIZ directly affects the 


time required to reach the optimum. 


B. PROBLEM FORMULATION 
The air-cooled heat exchanger is shown in Figure 10. A 


cross-flow arrangement with both fluids unmixed was chosen. 
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Cool air enters the heat exchanger at temperature ye ; 
Meessure p_, and constant specific heat c_ . The cool air 
makes one pass through the exchanger as feueiion over an isos- 
celes pitched bank of finned tubes. The air is heated by 
water, in single phase, at an entering temperature of Ty ; 


a 
and constant specific heat c_. 


W 
The analysis of the air-cooled, cross-flow heat exchanger 
centers about the first law of thermodynamics and on the heat 


transfer equation. These equations as they apply to the 


exchanger of Figure 10 are as follows: 


Q, =m Cour. -T  ) (12) 


Q, = mc «t. =f. } (p33) 


Q5 Be A OTL (14) 


wiere m=, T ana. mt. > “Lf are the fluid mass flow rates and 
a C5 W ha 

exit temperatures of air and water respectively. UL, 1s the 
true mean overall heat transfer coefficient based on the out- 
Side root tube area, A is the total heat transfer surface 
area of the exchanger used to compute Une and st 1s the 
mean temperature difference of the given exchanger. 

The object of the analysis, therefore, is to determine 


Q3, Qy and Qc. given an initial listing of values for the 


design parameters. The list includes the following: 
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Tubeside mass flow rate, lbm/hr 

Entrance temperature of hot stream, °F 

Exit temperature of hot stream, °F 

Beecific heat of hot fluid, BTU/lbm - °F 

Air mass flow rate, lbm/hr 

Entrance air temperature, °F 

Exit alr temperature, °F 

Entrance air pressure, psi 

Specific heat of air, BTU/lbm-°F 

Tube inside diameter (ID), in. 

Tube outside diameter (OD), in. 

fin height, in. 

mim thickness, in. 

Berm Spacing, in. 

Transverse Pitch, in. 

feng itudinal Pitch, in. 

Bank height, in. 

Bank width, in. 

Cross-flow arrangement 

Fin type 

Number of rows 

Number of passes 

Given heat transfer rate, BTU/hr 

Thermal Conductivity of tube material, BTU/hr.-ft.-°F 

Maermal Conductivity of fin material, BTU/hr.-ft.-°F 
Among the initial listing of design parameters above, 


there are parameters that are known and will remain constant 


a2 





macoughout the design problem. Also, there will be those 
parameters that are unknown and can vary, 1.e. design varia- 
bles. On the way to determining the various heat transfer 
rates, other information will have been computed. This inform- 
ation includes the objective function, constraining functions, 
and other design data, such as the number of tubes per verti- 
cal row. 

The optimizer will then manipulate the design variables 
in order to find an optimum, while at the same time, performing 


a heat balance, that is: 
where Q may be some given heat transfer rate. 


C. PERFORMANCE CALCULATION PROCEDURE 

With the temperatures, mass flow rates, and specific 
heats all specified in the listing of design parameters; 
whether they be constant or variable, the only unknown quan- 
tities on the right hand side of equations(12) through (14), 
are Oe A, and AT They will be determined as shown in 
Beoure 1. 

l. Mean Temperature Difference (MTD) 

For many flow arrangements, various approahces for 
determining MTD, mainly using diagrams, are available [36], 
which have proven very useful in manual design efforts. For 
computerized design, however, an explicit, approximate equation 


is desirable in order to achieve a fast, sufficiently accurate 
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calculation of the mean temperature difference of a given 


mow arrangement. 


Roetzel, et al. [37], presented such an approximate 
equation together with empirical coefficients for nine 
counter-current cross-flow arrangements as they apply to 
alr-cooled heat exchangers. 


Roetzel used the familiar equation for the MTD of the 


given flow arrangement AT: 
‘Vi Es AT (16) 
m 


where 4T, is the limiting case of pure countercurrent flow: 


m 
CE = ea i 
ae h c h Cc 
OTe = ile 2 Z i (17) 
(Th et. ) 
in 2 
(Ty pt ) 
2 i 


and F is a correction factor determined by a different set 
of coefficients for each flow arrangement. Roetzel reported 


m@emtOliowing function suitable for F: 


Mm n 


Mee tf - )}llCUl) a, (ly, )* sin(2i arctan R) (18) 
i=1 k=1 7° m 
where Va is the dimensionless LMTD. 
Mm 
AT, 
ey Te, 04 
iE 1 
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R = eee (20) 


and the coefficients ask of the approximating equation (18), 
were Calculated using a standard least squares estimation 
program [38], and are reported in reference [37]. The assump- 
tion that both streams were unmixed was used in their calcu- 
mation. 

When more than four tubeside passes are used, it is 
assumed that the heat exchanger has approached the limiting 
case of pure counterflow and F is set equal to one [2]. 

Having determined the MTD, the remainder of the analy- 
Sis procedure follows Roetzel's [39] general approximation 
method for determining the mean overall heat transfer coeffi- 
eient, Us for any flow arrangement while taking into account 
the variation of the heat transfer coefficients and the 
pressure drop with temperature and/or length of path. 

Before continuing with specific analysis procedures, 

a brief summary of Roetzel's general approximation method is 
me Order . 
The local overall heat transfer coefficient based on 


the outside root tube area can be written as follows: 





a 1 
eae A ents 77 :) (21) 
BGreg. Cee 3 wil 
A. H. 27KL Hn 
en Oat 


where Ne is the efficiency of the extended surface. 
In order to determine the individual convection heat 


transfer (film) coefficients, H and Hs according to the 
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conventional methods, the coefficients would be considered 
constant, and the necessary fluid properties for their calcu- 
lation would be evaluated at some mean bulk temperatures, 


aT ema TT. . 
C 


b b 
However, the film coefficients are not constant, but 


h 


vary with temperature and/or length of flow path. Roetzel 
has taken these variations into account with the use of 
corrected reference temperatures. Two sets of corrected 


reference temperatures are determined: Ty ; a and 
I ah 


dh a Therefore, for each set of corrected reference 


; é 
“a ie 
temperatures, the film coefficients are determined in the 
conventional manner using the reference temperatures in place 
of the bulk temperatures. 


With the film coefficients, Hi - A F H and H . 
I I ET aan 


two local overall heat transfer coefficients can be calcu- 
lated from equation (21), UL and Ul: 
Finally, the true mean overall heat transfer coeffi- 


cient is calculated as: 


ee (22) 


ea 
2. Determination of Reference Temperatures 
The reference temperatures for a pure counterflow 
heat exchanger must first be determined from [39]: 


= (Ty le ) 


: 1 
Hess = aa TPT =—7T: () T=) ye (TT =f 7) 
hy he Ay h, hy Co ho Cy 


(23) 
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2 i 
TT =e + (T =v 7) (24) 
h h h Cae (Tee) 
Le Z L 2 Ay Co h, 1 
AT. - (Ty -T ) 
Z ih 
a! =e foie =. on) (25) 
C Cc E c T -T )-(T -T_ +) 
I a 2 1 hy C hy Cy 
AT SIT (Ty tae ) 
1D = T + (T -T  ) ee meme 2s (26) 
e C S e 
AE i 2 il (T -T aa =a ) 
1 Z 2 iL 
where 
AT. = (Ty -T. ) + 18868 ; (oe -T y 41132 (olga 
Hi Z 2 1 
_ ay 
ATL, _ (T, -T. Za 2) ie -T 8868 (26b) 
1 2 Z Re 


Equations (26a) and (26b) were derived by Roetzel in references 
[40] and [41]. For the special case where the fluid heat 


Baeseaty rates, C, are equal, that is, 


m 
h Ph Cc Po 
the term 
Lilet. =T. ) 
J No Cy 
feel) —  ) 
No © Ny Cy 


of equations (23) through (26), with j3 = I and II, becomes 


-/8868 or .21132 respectively. 
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3. Correction of Reference Temperatures for Given 
Arrangement 


With inlet and outlet temperatures fixed, pure counter- 
flow yields the highest mean temperature difference. There- 
fore, for any other arrangement, the temperature difference 
would be smaller. Thus, the corrections are applied in the 


momLowiIng manner [39]: 


T = T' - y 
Ay Ay Ay 
T = T' -y 
Ary a Ary 
To = To tan 
ri I I 
o ae ee 
mb UT ba 


where Yh and ¥, are the temperature corrections of the 
k k 
hot and cold streams, respectively. The corrections are 


calculated as follows: 


ieee 7 AT 
m ios 


p =— AT. 
hi 5 : 27/2 
5 1 + (C,/C.) 
1- AT Ay 
w = AT 8 Sa aT Tad 
Cz 5 Se cr 
J 1+ (C./C,) 


The corrected reference temperatures are now used to 


determine the thermal conductivities and absolute viscosities of 


38 





fluids for later use in the calculation of the film coeffi- 
rents. thermal conductivity and viscosity data are usually 
presented in tabular form. However, for use on the computer, 
an explicit, simple approximate equation with temperature as 
the independent variable, was desirable. Water and air were 
chosen as two fluids that were likely to be involved in an 
alr-cooled heat exchanger design. Figures 1l, 12 and 13 
indicate that the thermal conductivities of air and water 
and the viscosity of air can be approximated by a second order 
polynomial. The viscosity of the hot tubeside fluid, water, 
must be treated specially, due to the following considerations. 
Calculations for the tubeside heat transfer coeffi- 
cient in the laminar, transition and turbulent regions all 


-14 [42]. In the past, an 


involve the correction, (u/u, 17? 
iterative procedure was required to determine the inside tube 
wall temperature in order to evaluate Mall’ Roetzel [43], 
has developed an iteration-free method for determining this 


correction. Roetzel's method assumes that the tubeside 


fluid's viscosity follows Andrade's Correlation, that is: 


Using the viscosity data of saturated water from reference 
[44], the coefficients were determined through regression 


analysis with the resulting equation: 


Zito. 7 104 


u = .01339 exp( 7 


te ) (27) 


So 





where temperature is in degrees Rankine and viscosity is in 
lbm/ft-hr. Equation (27) was plotted next to the experimental 
data with results that indicate that water does indeed satisfy 
Andrade's equation, see Figure 14. 
4, Uncorrected Tubeside Heat Transfer Coefficient 

The tubeside heat transfer coefficient is calcu- 
lated from one of three Nusselt-type empirical equations 
[42], as follows: 

For laminar flows, Reynolds Number < 2,100 
3 


ei 


L 


emanate (28) 


Ywall 


—— = 1.86 Re Pr ( 


For transition regions, 2,100 < Re < 10,000 





HD. DD: 
= = rRe*/3 ey (=e) 74) pri/3 (2 is 
u 
wall 
(29) 
FOr GUubobulent flow, Re > 10,000 
HD. 
== = 023 re? pri/3 (Hy +14 (30) 


WwW 


Therefore, before any calculations can even begin, 
the Reynolds Number, Re, must be computed to determine the 
[moe Or flow. 


The Reynolds Number will be calculated as: 





ee , mM 


A is supplied in the initial list of parameters and 


h 


Un is calculated from equation (27). This leaves only the 


cross-sectional flow area, AL to be determined, where: 


The number of tubes per pass, Bo is determined geometrically, 
having been given the tube bank height, number of rows, num- 
ber of passes, transverse pitch, diameters, and fin height 
initially. Figure 15 details the procedure of finding a 


Mathematically, 


n-ows re 
p = a = D = =z 1/P 


ie 


and would be a rational number. This rational number is 
then truncated to an integer number of tubes. 
With Re calculated, the Prandtl Number, 


€ 


1s computed from given and previously determined thermo-physical 
properties. 


The uncorrected tubeside heat transfer coefficients, 


femepetcnat is; H. , without the factor (u/u y 14 can be 
Calculated from the proper choice of equations (28) through 


(29). 


41 





5. Calculation of Wall and Associated Resistances 
Equation (21) can be written in a more general form 


as: 


1 
+ 
iL - oe 7 


where R., Rg and Ry are the inside wall and outside heat 
transfer resistances respectively. Additional resistances, 
such as contact or fouling, can also be added here. 


In comparing equations (21) and (31), the resistances 


Can be computed as: 





os a 
ee Ds He tee 
J a ee 
J 
eee knoe 7 Oe) 
Ae Sep (31b) 
wall 
_ ile 
| ae ae ee 
, @e.- 5 
J 
oe Calculation of Airside Heat Transfer Coefficients 


In order to have the optimization program play a 
Significant role in the selection of an optimized surface 
for a finned tube heat exchanger (which is one of the objec- 
tives of this project), an explicit equation for He involving 
tube and bank design parameters as independent variables is 
a necessity. 

In the past, comparison methods, as described by Shah 


Mao), were used to choose the best surface from a list of 
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surfaces for which experimental heat transfer and friction 
data existed. The data is presented in graphical form, 
where Friction and Colburn Factors are plotted as a function 
of Reynolds Number. With J, mn, a and Pr known, the film 
coefficient can be computed. 

Therefore, in previous heat exchanger optimization 
programs, a given surface described by its pitch, outside 
fin diameter, fin thickness, fin spacing, and outside tube 
diameter, had to be chosen beforehand. After the surface 
configuration had been selected, expressions for f and J 
were obtained by fitting polynomials to the experimental 
data described earlier. 

Briggs and Young [47] provide a means of getting 
past this obstacle with an improved convection heat transfer 
correlation for air flowing across triangular pitch banks of 
high finned tubes. 

Briggs and Young expanded the work of Ward and Young 
[48], to cover a total of 18 differently configured finned 
tube banks in order to determine the effect of fin thickness 
and tube pitch on the airside heat transfer coefficient. The 
heat transfer data for the high-fin tube banks were correlated 
Eengive : 


x O 
Nu = —_— = (1378 re’ (+8 p,1/3 (3) 7478 (32) 


where s is the distance between adjacent fins and 2 is the fin 


height. Equation (32) is based on tubes having a wide range 
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of fin heights, fin thicknesses, fin spacing and outside 
tube diameter and can be used to predict A for a bank of 
tubes six rows deep. Figure 1l of reference [48] is used to 
Serrect HS. for banks of other than six rows. 
os erieai ation of Fin and Surface Efficiencies 

In order to calculate the outside heat transfer 
resistance, which will be used to calculate the correction 
for the tubeside heat transfer coefficient and finally, the 
local overall heat transfer coefficients, the extended sur- 
face efficiency, Nes must be computed. 

The surface efficiency accounts for the temperature 
@rop from the root to the tip of the fin, due to the thermal 
resistance of the fin material. Thus, even though the heat 
transfer has been increased by the additional area of the 
extended surface, the area is not as effective as if it were 
to be all at the root temperature. 


The surface efficiency can be expressed as [49]: 


As 
Ne — eta (> e) 


where ¢, the fin efficiency for a radial fin, can be found 


from [49]: 


2 Lo ee a! ~ ky (mr) I, (mr) 
2 Zee ite eK 
=, ) OO Oo 


J 
: ; (mr ¢) + I, (mr) k_ (mr) 


i: 


The finned area, Ags and the total heat transfer area, A, are 


computed as follows: 





With Ro calculated from equation (3lc), the correction to 
the tubeside film coefficient can be made. 

8. Correction of the Tubeside Heat Transfer Coefficient 

With the tubeside heat transfer coefficient dependent 

on the wall temperature, the dependence has either been 
neglected, or the wall temperature has been calculated with 
an iterative technigue in the past. Roetzel [43], has pro- 
posed an improved iterative-free method for finding the "Seider- 


14 


Tate" correction, (u/u )°"", when the tubeside fluid 


wall 


obeys Andrade's viscosity correlation. 


From Roetzel's formulation: 








2 
Ma ke 
(my FF 2 - De (pt vi? 
wall 5 J 
where 
* 
“wall, 1 i= T. 7th. 
u. = .007 in( a ee 7] 
J Pe Ba | Be are 
J J 
* 
cay lea 1 ae cae on 
eg ee) aa s~ na = 1 Booare 
J J j wall J 45 


45 





where 


—_ ! 
25 “i a Sl op or 


4 


(Note that all temperatures are in degrees absolute.) 

All the parameters on the right-hand side of equation 
(21) are ROA Aiable, Therefore, the two reference over- 
all heat transfer coefficients, Us and Usy: can be calculated. 
The mean overall heat transfer coefficient follows easily 


from equation (22). 


With Ua! Q. can be calculated, with the heat balance 

to be performed by the optimizer. 
9. Calculation of Air and Tubeside Pressure Drops 

The final calculations before computing objective and 
@emocchaint functions for the optimization problem involve 
the pressure drops in the heat exchanger. 

The basic equations that will be used for the calcu- 
lation of the pressure losses are as follows: 


a. Tubeside [42] 


2 2 
f. m, wn (n_-1) mm, 
Ap, = Sg ee P (33) 
2g_A o0 D. (—) CoC. p 
C Xm oe Cc = 
where z = .14 below Re = 2100 and z = .25 for Re greater than 
ZOO. 


b. Airside [50] 
fs n. m_* 

‘DP. = We (34) 
fe te Jo P 


where Age is the minimum flow area. 
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The friction factor for the tubeside flow, fi, is 
taken from Figure 9.5 of reference [42], assuming fully 
developed flow. For use on the computer, an explicit expression 
mor f. was obtained by fitting a line and an exponential to 


the experimental data of Figure 9.5. This follows: 


Re < 1000 f£, = .5/Re 
Re > 1000 cee oor 2: Re °©674 


Just as in the case of the airside film coefficient, 
for surface optimization on the computer, an explicit equation 
Berethe airside friction factor, foe is a necessity. Robinson 
and Briggs [50], presented such an expression for fy for air 
flowing across triangular pitch banks of finned tubes. Robinson 


and Briggs' work closely parallels that of Briggs and Young 


[47]. The Robinson-Briggs Correlation: 
I Pe 
f - 18.93 Re OkO (_t) SPT! (t) -515 (34a) 
O D Pe 


covers the range of tube sizes and pitches used in air-cooled 
heat exchangers [50]. 


Therefore, with Re. and Re the four reference 


i ee 


Pressure drops, Ap; , and AR, , May be computed from equations 
: : 
(33) and (34). Following Roetzel's general approximation 


method: 
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For a gas, an additional correction is needed because the 
density in equation (33) is strongly dependent on pressure, 
which is changing through the exchanger. Using the inlet 


pressure as reference: 





Ap Ap 
—— “I, il tes ae, 
© Uy Urt ae 
Zap. 
_ _ Ol, 2 
AP, = Py EL (el Py ) ] (35) 


All the necessary information from an analysis view- 
point has now been calculated. Functions needed for the 
numerical optimization process shall follow. 

ie Objective and Constraint Functions 
The objective functions available for minimization 


are defined as follows: 


a. Volume = wh[D, ~ (n -1)P, cos 0] 
where 
© = arcsin (P,/2P,) (36) 
y 2 
(p -~D “)r 

wo f O me De 
b. Heat Transfer Area = i Sse + (lL 5 O 

Ap. += m™ 


c. Air Horsepower = 


48 





ole 


e. 


Airside Pressure Drop 
Tubeside Pressure Drop 


Limitations were imposed on the following inequality 


constraints in order to keep the design within practical 


physical bounds: 


Ae 


The diameter ratio, 


De 
DRATIO = = 
D 
O 
must be kept reasonable. This can also be accomplished 


to some extent by placing side constraints on the 
design variables, 2% and Da: See Figure 16. 


The optimizer must be prevented from driving the 
tube thickness, 


TUBTH = ee oye 


me Zero. 


The tubes must be kept from touching in both the 
longitudinal and transverse directions, 


LOUGH NI — 2 ai= 


FOUCHT 


I! 
O 
I 
tg 


TOUCHL and TOUCHN, therefore, must be Kept below zero. 
Reasonable temperature profiles must be maintained 


at both ends of the exchanger. See Figure 10-7 of 
Ref. [44], 


AR Olg sl 


it 
FH 
: 
KF 


Prone 


tl 
a 

] 
FF 


that is, PROFH and PROFC, must be negative. 
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The number of tubes per vertical row, VROWR, shall 
not be allowed to go below 2. 


The free face area, that is, the minimum flow area 
for air, muSt obviously be greater than zero, 


DELSFF = [projected tube area] - hw 


where the projected tube area, STOTAL, is 


eo 





STOTAL = N_[D_w + 
ro 
and N. is the number of tubes per vertical row. DELSFF 
must Be less than zero. 


The airside and tubeside pressure drops must be kept 
within design constraints. 


From experience, the argument of the square root of 


eagtation (35), 


0 
ARGS =) i= 24p' 
Py 


has been driven below zero. It must therefore be 
constrained. 


TO Maintain an isoscele pitch bank, the angle, 0, 


as shown in Figure 17 and as defined in eq. (36), must 
be constrained. One such constraining value is: 


oe = arcos (D-/2P, ) (37) 


The arguments of the arcsin and arcos of eq. (36) and 
(37) must be constrained from going beyond l. 


The heat balance described by eq. (15) is performed 
by three equality constraints: 


QRATIO = Q./Q (38) 


QRATIOL = Q,/0, (39) 
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QORATIO2 = 04/2, (40) 


All constraints are set equal to l. 
Besides the constraints described in paragraphs a. through 
j., above, lower and upper bounds are placed on the design 


variables to assure a reasonable design. 


Sek 





iV URESULIS 


A. BACKGROUND 

Case studies were chosen as the best way to test the 
Seapabilities of the program for Heat Exchanger Design using 
Numerical Optimization (HEDSUP). The design problems posed 
were made as realistic as possible. 

1. Capabilities 

HEDSUP currently has the capability to design for 

Nine different configurations of triangular pitch banks of 


finned tubes: 


Dini ROW, 1 PASS 
Pete 2 ROW, 1 PASS 
ia oS ROW je AS > 
TYPE 4 - 4 ROW, 1 PASS 
ibe 2 RON, 2 PASS 
IVene owe 96 ROW, 3 PASS 
twit? 9 94 ROW, 2 PASS 
TYPE 8 - 4 ROW, 2 PASS 
ieee PURE “COUNTERE LOW 


fieeeetO will include exchangeérswith a configuration of n rows, 
meeasses, where n can go from five to 20. 

The banks must be constructed of high-finned tubes 
(2 > .0625 in. [42]) with the fins having a rectangular pro- 
file of constant thickness. Additional profiles can be inserted 


into HEDSUP quite simply, provided that its fin efficiency can 


a2 





be expressed explicitly as: 
n =e Cl een, 1) 


See Subroutine FINEFF of the program listing, Appendix D. 

At present, HEDSUP can provide the design parameters 
for an air-cooled heat exchanger optimized for any one of the 
following design objectives: 

(a.) Minimum Volume 

(b.) Minimum Heat Transfer Surface Area 

(c.) Minimum Air Horsepower 

(d.) Minimum Airside Pressure Drop 

(e.) Minimum Tubeside Pressure Drop 
Additional design objectives can be used, provided that they 
Can be expressed explicitly as a function of the design varia- 
bles and they are added to the common block. It should also 
be pointed out that any deSign variable may simultaneously be 
a design objective as long as it conforms to the restrictions 
Of both. For example, an exchanger may be designed for mini- 
mum bank height. 

The airside fluid is restricted to dry air. The tube- 
side fluid is presently limited to water in single phase. 
Other tubeside fluids can be used by HEDSUP, provided that 
their visocisites obey Andrade's Law and the fluid's thermal 
conductivities can be expressed explicitly as a function of 
temperature. The fluid's specific gravity would also have to 


be placed in the denominator of eq. (33). 
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Ee CASE STUDIES 
fm Case One 
a. Problem Formulation 
An air-cooled heat exchanger is to be designed 
for minimum volume with a heat transfer rate of 10,000,000 
BTU's per hour. Water is to be cooled from 200°F to 125°F. 
Dry air will enter the exchanger at 95°F and leave at 130°F. 
Specifications call for a fan that can produce a pressure 
difference of two inches of water. 
b. Design Variable Framework 
From the list of design parameters in Section 
III.B, the design variables must be singled out, identified 
to COPES, and given side constraints. All parameters must 
be given an initial value. Only the values of the design 
variables will change. 


Assuming constant specific heats, 


Cc eee. © Bre on! 


Cc = .24 BTU/1lbm- °F 


the mass flow rates of both fluids can be determined from 
eqs. (12) and (13), because the heat transfer rate and tem- 


perature differences are given, 


ar: ; 
m= a 133333. ibm/hr 
Ph 
a eo oavicne aon /he 
a a AT. c 
C 
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A cross-flow arrangement, fin profile, tube 
material, and fin material must be chosen. 


The design variables for this example are, there- 


mole : 
.232 < D, < 2.325 in. D,* = 2.0 inches 
~ aE D eS alas D.* ee Caes 
POGOZ5 tm. < {<2 gi = .46 inches 
01 <¢€ < .0235 in. t+ = .023 inches 
MOM Seem a 25 IN. st = .1l11 inches 
0.0 <P. < 4.0 in. Pp = 2.125 inches 
00° 22< Py Sean eT P,* = 4.0 inches 
0.0 <w < 500 in. wt = 490 inches 
OOo fh < 500 “in. ha = 350 inches 


The side constraints on the design variables are of a practi- 
cal nature with the exception of the lower bounds on fin 
height. Recall that the use of eq. (32) is restricted to 
high-fins. High-fins will also tend to keep the fluid unmixed, 
which was an assumption used when defining the coefficients 
ae kk! used in eq. (18). 

c. Constraint Framework 


From the problem statement, the airside pressure 


drop must be less than two inches of water or .0722 psi, 


a 





Oe < dp, < SOM 2 OS 1 

From a practical standpoint: 
One Cale]. 3 
ia <—seDRATLO 9 92.5 
LOL¢CueaeTUBIH <9, 28 “anches 
-© aS ~14 psi 


The equality constraint 


where 
Q,. = UA AT 
Oue— elo ouo- O00 — Btu/hr 


will satisfy the heat balance. 

The bounds on TOUCHN, TOUCHL ARG5, ARG7, ARG8, 
DELSFF and VROWR were discussed in Section III.C.10. Con- 
straints on PROFH and PROFC were unnecessary because all 
temperatures were specified in the problem statement. 

As with design variables, constraints must be 
identified to COPES by location in the common block. See 


the global catalog, Appendix A. 
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ad. Methodology 

Ideally, a three dimensional design matrix can 
now be constructed of optimum exchanger designs with minimum 
volumes. The matrix would be constructed by first holding 
the tube and fin materials constant and varying the configura- 
tion, 1.e. Type l, Type 2, Type 3, etc. Next the tube material 
would be varied with the fin material and SRGAeee configura- 
tion held constant and so forth. However, for case study one, 
the tube material will be chosen as copper, k = 200 BTU/ft-hr-°F, 
and the fin material will be aluminum, k = 118 BTU/ft-hr-°F. 

Also, in order to simulate an actual trade off 
study, the constraint framework will be fixed throughout the 
individual case studies. 

Problems arise in constructing the matrix when 
trying to determine the true minimum volume design for each 
configuration. Unfortunately, the choice of initial design 
parameters (starting point), coupled with the input parameters 
for ALMM, will sometimes lead to entirely different optimum 
designs with volumes differing by over 100%. 

The ALMM parameters include the initial multiplier, 
ee, ene multiplication factor, CMULT, and the maximum multi- 


plier value, CCMAX. Experience has shown that setting 


Chis = 2.0 


CCMAX 1000. 


will suffice for almost all problems. However, there is much 


"artwork" involved with the choice of CC. From experience, 
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an initial multiplier of 10 works well when starting far 


from the equality constraint, i1.e.: 
Tome SF ORATEO < (1.2 


However, when approaching the heat balance, i.e., the equality 
constraint, a CC = 100 is necessary for the heat balance to 
converge. Equation (15) is considered satisfied when the heat 
transfer ratios (eqs. (38) through (40)) are less than 0.1 
percent. 
Therefore, it is obvious why the chosen initial 
design is so critical. Together with the choice of CC it 
will determine how and to what design the optimizer will con- 
verge. As an example, see Table 1. Notice the calculated 
heat transfer for the initial design. This value is the product 
oi Ut A, and AT. calculated by ANALIZ using the initial design 
parameters, some of which are mere estimates. Recall it is 
the job of the optimizer to vary the design variables in order 
=O bring oF equal to QO, Q, and on and atthe same time minimize 
the objective function. 
In order to remove some of the "artwork" and 
try to ensure a true optimum design, the following procedure 
is recommended to generate the design matrix: 
(1) Begin with a Type 2 configuration; input the initial 
design and constraint values enumerated in sections 
(a) through (d) above; let CC = 10.; execute. 
(2) If the heat balance of the resulting design has not 
converged, but is within 20%, use the design results 


as a new Starting point and set CC = 100 (if the heat 
balance is not with 20%, let CC = 10.). 
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TYPE 7 EXCHANGER CONFIGURATION OPTIMUM DESIGN 
























CETIMUM 
DESIGN 


OPTIMUM 
DESIGN 


PE As 
DESIGN 


INITIAL 
DESIGN 


DESIGN 
VARIABLES 
















inches (Plas 
ches pees 
inches .294 
inches HOS 
inches .08 
inches Looe. 
inches ieee 99 
inches eo 
inches 253% 






fee, BTU/hr oS O24 ino foo, 225. 120,600,352. | 9,999,ell. 


velume , fe 23s i2 P22 BO S.22 i520 5 





Table l 
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(3) Repeat step (2) until convergence. NOTE: CC may 
be adjusted up to 150 when approaching convergence. 
The design should converge following the use of 
cC = 100. If too much adjustment of CC is necessary, 
reaching the optimum from the starting point is 
unlikely. 


(4) Ensure a minimum design by beginning step (1) witha 
different initial design. 


(5) After finding an optimum design for a Type 2 con- 
figuration, use that design for the starting point 


for a Type 3 configuration. This assures a reason- 
able starting point, probably close to the optimum 
Bor Type 3. 


(6) Repeat steps (1) through (3) for a Type 3 configuration. 
(7) Repeat steps (5) and (6) for the remaining configuration. 
e. Design Matrix 
Hieseaesign Matrix 1S presented in Table 2. The 
optimum design is a Type Four configuration and is shown in 
Figure 18. Typically, for this case study, when starting far 
from the final design, COPES/CONMIN would require approxi- 
mately 1700 calls to ANALIZ to reach an optimum. However 
when beginning from a reasonable starting point with CC = 100, 
only 600 calls were needed. Note that each call to ANALIZ 
requires approximately .06 seconds of CPU time on an IBM 360/67. 
Z. Case Study Two 
a. Problem Formulation 
An air heater is to be designed to fit into a 
space 8' X 24' X 4', The heat exchanger is to heat 1,000,000 
lbm/hr of dry air from 75°F to 130°F. 256,000 lbm/hr of water 
at 200°F is available. Design the heater so that the required 


air horsepower is at a minimum. 
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b. Design Variable Framework 


Assuming constant specific heats, 


C = 1.0 BTU/ lbm=°F 


e = .24 BTU/lLbm- °F 


eqs. (12) and (13) will yield the required heat transfer 


rate and the outlet water temperature. 


DC Glentreee =r) = is.2<10° BTU/hr 
Ps -5 al 
my “p, "hy en @) 
ee «OP 
2 my, 
Therefore, the design variables are as foliows: 
eee eS) 2-325 In. Bae SO NO: aa 
25 < D <«< 2.5 in. D7* = 2.5 in. 
© © 
.0625 in. < 2 < = g~* = .46 in. 
Ol < t  < .0235 in. t* = .023 in. 
Mise SS = «6125 in. s*> = .111 in. 
i 
0.0 < Py < © PS = ile 5 heats 
i, 
OO < Py < Py =— G00 2h. 
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c. Constraint Framework 
When designing for minimum horsepower, the opti- 
mizer will naturally try to drive the design to a maximum 
volume in order to reduce airside pressure losses. Therefore, 


it is reasonable to constrain volume as follows: 
O < Volume < 768 oo 
The other constraints are: 
O20mr< 0 <eel. s 
iO ome Oka eo: < 82.5 
Ses TUBTH < .18 ain. 
QO < 4p, -14 psi 


and for the heat balance: 


e 


2-1 


Q 


d. Design Matrix 

With the tube and fin materials fixed, as in case 
one, the matrix iS presented in Table 3. The optimum design 
is a Type 4 configuration, as shown in Figure 19. In case 
study two when starting far from the optimum, the optimizer 
called ANALIZ approximately 1900 times. When starting close 
to the final design, for example using the design for a Type 3 
configuration as a starting point for the Type 4 design, 


COPES/CONMIN only required 596 calls. 
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Wee  CONCEUSTONS 


The intent of this investigation was to couple an analysis 
program with a numerical optimization scheme, COPES/CONMIN, 
to produce a complete, detailed design program for an air- 
cooled heat exchanger, (HEDSUP). In addition, the analysis 
program was to be written such that: 


1. the variation in the film coefficients with temperature/ 
length of flow path would be taken into account 


2. the surface would be optimized 


3. it would be iterative free and thus minimize the CPU 
time required during an actual trade off study. 


The results from test cases using ANALIZ coupled with 
COPES/CONMIN in its present form were unsatisfactory. Although 
COPES/CONMIN could optimize the objective function satisfying 
the inequality constraints, a reliable heat balance could not 
be obtained. The solution to this problem was the addition 
of the ALMM option to COPES. In this way, the method of 
feasible directions, which works best with inequality con- 
straints, was used to satisfy the inequality constraints. 

The multiplier method, which works best with equality con- 
Straints, was used to perform the heat balance. 

The results of the case studies show that HEDSUP will 
yield reliable designs for various design objectives and 
problems with only some trial and error application of the 
initial Lagrange multiplier. However, precautions must be 


taken to overcome the relative minima that plague this design 
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problem. Table 1 shows vividly the problem of relative 
minima with one "optimum" design having a volume over 100% 
greater than the "true" optimum. 

The value of numerical optimization in a deSign problem 
of this size cannot be overemphasized. For example, in case 
study 1, the problem is taking place in a nine-dimensional 
design space and intuition on how an optimized design "should" 
turn out is quickly lost. Figure 20 helps to illustrate this 
point. When only varying two design variables, h and w, 
Figure 20 shows that, when beginning from the initial starting 
point with a Type 7 configuration, a design satisfying both 
the equality constraint and the air pressure drop constraint 
could never be found. However, Table 2 indicates how the 
numerical optimization routine has varied the other seven 
design variables in order to shift the constraints and yield 
an optimized design. 

The results show that AHDOP did vary the surface design 
Meaerables: D., D, &, t, S, P. and P 


aL Oo t ie 


an optimum heat exchanger. This capability of surface opti- 


in order to produce 


mization is dependent upon the use of the Briggs-Young and 
Robinson-Briggs Correlations, eqs. (32) and (34a), respec- 
tively. The reliability of the correlations as compared to 

the "conventional" method is questionable. Actual experimental 
data for a particular tube and pitch will always be the most 
useful in predicting pressure drop and film coefficients of 


across banks of finned tubes. However, the correlations 
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~~ 


mentioned above cover the ranges and pitches used in air- 
cooled heat exchangers, and should therefore be sufficiently 


accurate in predicting He and AD .° 
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VI. RECOMMENDATIONS 


In addition to the insight that this investigation has 
given into the generation of an automated air-cooled heat 
exchanger design, it has also generated an awareness of this 
investigation's shortcomings. As mentioned in the review of 
previous work in this area, each of the optimization methods 
has its own limitations; none is completely general. Pre- 
sented herein are recommendations for improving upon and 
furthering development of HEDSUP. 

1. HEDSUP should be expanded to include the capability 
for two-phase tubeside fluids. Mott, et al. [2] discusses a 
method involving two-phase tubeside fluids that would be 
compatible with HEDSUP. The modular design of ANALIZ will 
eelaethnis effort. 

2. Research in the area of numerical optimization using 
discrete variables would benefit HEDSUP immensely. With dis- 
crete variables, the design of the exchanger could be accom- 
plished with "off-the-shelf" materials. At present, the use 
of the optimizer is restricted to continuous variables with 
continuous first derivatives. The ability to work with dis- 
crete variables would also eliminate the need for the design 
matrix. The optimizer could optimize for type configuration, 
fin profile, fin material and tube material. 

3. Additional research with ALMM is needed to remove the 


"artwork" involved with choosing an initial multiplier and 
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thus increase reliability and hence reduce CPU time. The 
research should be concentrated in two areas: 1) debugging 
the new optimizer with AFMM and 2) scaling of design variables. 
4, The addition of cost as an objective function would 
increase the attractiveness of HEDSUP. Mott, et al. [2], 
and Fontein and Wassink, [18], have presented much useful 
marOormation in this regard. 
5. Mechanical constraints such as tube bursting stress 


and tube vibrations should be included in HEDSUP. 
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Calculate Log Nean 
Temperature Difference 


Calculate Cerrected 
Reference Temperatures 


. Calculate Number of 
Tubes & Tubeside Flow Area 


Compute Uncorrected Tubeside 
Heat Transfer Coefficient 


Calculate wall & Associated 
Resistances 


Compute Uncorrected Airside 
Heat Transfer Coefficient 


Correct Airside Heat 
Transfer Coefficient 


Calculate Fin & Surface 
Sfficiencies 


Calculate Fin & Total 
Heat Transfer Areas 


: 


Calculate Airside 
Resistance 


Correct Tubeside 
Heat Transfer Coefficient 


Compute Rererence Mean Over- 
all Heat Transfer Coefficient 


: 


Calculate True Mean Over- 
all Heat Transfer Coefficient 


Calculate Air & Tubeside 
Pressure Drops 


eee Constraints4 
Creccive sunctien 


Figure l 
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APPENDIX A 


HEAT EXCHANGER DESIGN USING NUMERICAL OPTIMIZATION (HEDSUP) 
USERS MANUAL 


In order to execute HEDSUP, it is necessary to provide 
formatted data for COPES [25], followed by formatted data 
for the ANALIZ portion of the program. This section defines 
the data which must be supplied. The data is submitted in 
segmented blocks. All formats are alphanumeric for TITLE 
micmeND cards, F1O for real, and I10 for integer data when 
supplying COPES. For ANALIZ, formats are F14 for real and I10 
for integer data. 

While the COPES user's sheets define the data in formatted 
fields of ten, the COPES manual does provide means of simpli- 
fying this task through unformatted data input. 

The included GLOBAL CATALOG defines objective functions, 
design variables and constraints along with their location 


for ease in compiling the necessary COPES data. 
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14 BANKW Bank Width, in 


Bank Height, in 
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[39 | Pek _|Aireide Pressure prop, pei 
[20 [Pew _|tubeside Pressure Drop, psi 
wi 
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- 4 BTU/hr 
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i. TUBTH Tube Wall Thickness, in 
eee TOUCHN Tip-to-Tip Clearance, Transverse 
TOUCHL Tip-to-Tip Clearance, Longitudinal 
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